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Abstract 

The weak field average magnetic susceptibility of square shaped mesoscopic 
conductors is studied within a semiclassical framework. Long semiclassical 
trajectories are strongly affected by static disorder and differ sharply from 
those of clean systems. They give rise to a large linear paramagnetic suscep- 
tibility which is disorder independent and in agreement with recent experi- 
ments. The crossover field to a nonlinear susceptibility is discussed. 
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Recent develop ements in the theory of thermodynamics on the mesoscopic scale have 
brought to the focus of attention the question of the average magnetic susceptibility of an 
isolated mesoscopic counductor, (x(H)). It has been shown that as the number of electrons 
in such a system is independent of the value of the applied magnetic field, one has to 
evaluate this quantity within the canonical ensemble. Upon energy averaging (or averaging 
over system size), one may discard oscillatory grandcanonical contributions, obtaining 

(x{H)) = -^^ f2 {8N\H))\-, (1) 

where (5N 2 ) is the typical (field dependent) sample to sample fluctuation in the number of 
levels below some effective chemical potential ~p IJ. Efforts have been concentrated on the 
calculation of (x) in the diffusive regime, but have been extended as well to the limit of 
clean systems (where the only scattering is due to the sample's boundaries) [0-0]. With the 
regime 

L > L (2) 

made accessible experimentally ||, the belief is that such theoretical studies may be em- 
ployed to interpret the experimental data. Here l e \ is the elastic mean free path, while L is 
the system size. 

From another point of view, it has become clear that the association of all systems sat- 
isfying Eq. (H) (perfectly clean or containing disorder) with a ubiquitous ballistic regime is 
too naive []9|JT0[1 . Let us consider a finite disordered system of an integrable geometry (e. g. a 
square). For disorder weaker than that which gives rise to diffusive behaviour, one has to 
compare not only L with l e i, but also the average level spacing, A, with the inverse elastic 
mean free time h/r e i. We may now distinguish between two weak disorder regimes satisfying 
Eq. (0): the ballistic regime ( h/r e i > A) and the clean (or perturbative) regime ( h/r e i < A 
or, equivalently, l e i > L(kpL) 1 . This regime is, presently, experimentally inaccessible.). 
For the latter, as disorder constitutes only a small pertubation on top of the clean spectrum, 
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one may argue that thermodynamic quantities are practically those of a perfectly clean sys- 
tem (of the same geometry). This is, however, not the case for the former, ballistic regime. 
It has been shown recently that elastic mixing of the levels give rise to non-trivial level 
correlations JIIJ (at H = 0). It is therefore natural to ask how thermodynamic quantities, 
notably the magnetic susceptibility, behave in this regime. 



In the present work we address the problem of the average weak field susceptibility 
of (integrable) disordered ballistic quantum dots. Employing a semiclassical picture [12| 



we find contributions to (x{H)) which arise from long trajectories that are sensitive to 
elastic scattering. However, due to subtle cancellations that occur within this framework, 



contributions to the zero field susceptibility do not depend on disorder p3| . We find that 
this susceptibility is given by 

X (H = 0) = +\ XL \ak F L . (3) 

It is paramagnetic and includes an enhancement factor akpL with respect to the Landau 
susceptibility, where the numerical factor a is estimated below. Our results are in qualitative 
agreement with the experimental data of Ref . |Sj . Our approach is closely related to the de 
Gennes Tinkham "method of trajectories", and some of our results bear close formal sim- 
ilarity to a method, discussed by Beenakker and van Houten, of treating weak localization 



corrections in restricted geometries [11,13 . 



Consider a possible multiply reflected semiclassical trajectory within our square geome- 
try. The same trajectory may be described within an extended zone scheme (Fig. la), where 
boundary scattering does not occur. We next consider a uniform perpendicular magnetic 
field in the z-direction. Within the Landau gauge the vector potential is A = —Hy~k. The 
vector potential in the extended zone scheme is shown schematically in Fig. lb. According 
to this scheme, the electron is moving in a staggered field |TJ| , which may be written as sum 
of its Fourier components 



16FLx ~ (-l) n „ . Tree . 
5 — Y 7 —tt 1 ^7 sin 2m +1 — ) 

sin((2n + l)^). (4) 



Imagine now a random walk trajectory (in the extended zone scheme) of total length 
C, starting from point (x ,y ), and consisting of C/l e i uncorrelated segments of length l d 
each. The end point of this trajectory may be mapped back onto the reduced zone scheme. 
For long trajectories (C > l e {) we assume that ergodicity holds and that the image of the 
end point in the reduced zone scheme is uniformly distributed. In the absence of an applied 



field, the probability of returning to the origin at time t is given by (see e. g. [|T^]) 

where the factor 2 in front represents an enhancement due to an interference of a trajectory 
with its time reversed image (a "Cooperon contribution" ). The two-level spectral correlation 
function, Yz(AE), is given then by the fourier transform of fll2f , and finally the level 

fluctuation is given by 



(5N 2 (AE)) = 2 / (AE - s)Y 2 (s) ds 
Jo 



A , K 2 p(t), ,AEt,, 

Our goal is to find the field dependence of pit), hence (SN 2 ), then employ Eq. ([!]). We first 
note that considering long trajectories, i. e. time scales larger than the elastic mean free 
time T e i, allows us to evaluate reliably energy windows such that AE < This, in turn, 
may be employed to evaluate the weak field susceptibility. (In other words, long trajectories 
affect the small energy contribution of the level-level correlation function to (5N 2 ) [[13| , 16f| ) . 
Estimates of the field scales for which the present treatment is valid are given below. Upon 
application of a magnetic field there is an extra phase that may be associated with the 
amplitude of the jth returning trajectory, (f)j = ^ A • c?x, where the integral is taken 
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along the jth trajectory. To calculate p(t) one has to take the square modulus of the sum 
over all trajectories corresponding to time t. Upon disorder averaging, only terms that may 
be represented as a product of an amplitude with its complex conjugate or with the complex 
conjugate of its reversed survive |I2| , ^7| . The factor 2 (Eq. ([)])) is now to be replaced by 
a factor (1 + cos 20,) to be averaged over all trajectories j . We note that following 
the disorder averaging, the dependence on random (trajectory specific) phases disappears. 
Our remaining task is to evaluate the distribution function of 2<fij. We first assume that the 
magnetic phases (contributions from line integrals) accumulated at different steps of size 
hi of a trajectory are independent and identically distributed, with a distribution function 
Pi el ((j)). The average over cos20j is then given by |H| 

C/l e l 



Re 
Re 
Re 



-oo 

1 r°° 

2lX J— OO J —TT 

— 7T 2tt 



1 £/lel 

(5(0 - 0(0))e* 2 * d9)d(<j>) 



(7) 



where 9, the scattering angle at a given step (cf. Fig. la), is assumed to be uniformly 
distributed (a consequence of isotropic scattering). The function 20(6*) (associated with 
a segment of length l e i) is evaluated by inserting Eq. (g) into the expression for the line 
integral. Approximating A x (x, y) by the m = n = terms in Eq. @), we obtain 

IHL r 1 ^ 



-2 ■ 2tt- 



1X- j <Pq JO 



ds cos 9 



(7T 7TS \ 

-^{xq - y ) + — (cos6» - sin #) J 

(7T ITS 
— (x + yo) + —(cos6 + sin( 



(8) 

where (xo,yo), (xo + hi cos 9, yo + l e i sin 9) are the endpoints of the segment at hand. We 
next write 



p(t) = -p (l + ReC(#) £/ ' ei ; 



(9) 



Accounting only for the field sensitive (Cooperon like) contribution Ap(t), writing £ = vpt, 
hi = vpT and introducing a cutoff factor [|l9j, 7, we obtain 



Ap(t) = ~PoC^e-£ . (10) 
Substituting this into Eq. (||), we obtain that the flux sensitive part of SN 2 (AE) is given by 



(AE) 2 + (Mn 1 7 + 1 ) 2 
-lni + 7) 2 

1 Tel C ' ' 



5iV 2 (A£) = p h 2 In ' Vr f ; _ C N9 (11) 



To evaluate the very weak field behaviour we may expand 

fl^.l./^^El-^. (12) 



Averaging over 2o,2/o we obtain 

c = — ds da' j d9 cos(tt\/2(s / - s) sin#) , (13) 

J-lel/L J-hl/L J-7T 



7T 5 



which, following a few straightforward steps, yields 



512 riei/L , , , r- .X 



c 



du J {\/2im){-f - u) , (14) 



7T 4 Jo ' " UV ' " M L 

where Jo is a Bessel function. As ^ 3> 1, this integral may be approximated by 

256^2 Z e * 



c 



(15) 



7T 5 L 

The most remarkable feature of Eq. is the proportionality of c to Z e /. The main contri- 
bution to the quadratic correction cH 2 comes from segments oriented at angles 9 ~ ±7r/4, 
for which (0 2 ) oc / 2 ; . However, the angular width around 9 = ±7r/4 corresponding to 
such exceedingly large contributions is ~ y-. rendering c ~ l e \. It follows then that 

'•el 

£ p tt2 t 4 

n e i ~ 1 _ -^c T3 , rendering the quadratic term, hence the zero field susceptibility, 
disorder independent. Employing Eqs. (HD and fllTf) , and considering an energy window 



AE ~ — 3> A, we obtain 

Tel 



A +3 2cL 4 1536^ A 



( X (H = 0)> = -po^ — -2 = r \XL\k F L 

7 ^ 00 ^" 7 

~ 0.23-|xl|W, (16) 
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where for spinless electrons xl = — of^i - F° r \ s P m electrons %■> as we ^ as \Xl\-, should 



be multiplied by 2 |23| . The range for which this very weak field analysis may be employed 
(i.e. the regime of linear susceptibility, eq. (0)) is found by (i) examining the range of 
validity of Eq. ([12]) , i.e., by requiring that ( may be expanded to quadratic order in H 2 and 
(ii) by requiring that the field dependent correction in Eq. flTID be smaller than 7! These 
conditions read 

H < (flux in the sample < ±<j) ) [24] (i) 



Typically, inequality (|l7|ii) is stricter, and yields the dependence of the crossover field on 
both the elastic and inelastic scattering. Beyond this range the susceptibility is not constant, 
and may even change sign (cf. Fig. 2). We next recall that our semiclassical considerations 
apply to energy intervals (close to 6f) smaller than We now find the largest field for 
which this picture may be employed to evaluate the total if -dependence of (5N 2 ) (hence 
x)- The condition is that ((2(f)(6)) 2 ) < 1, yielding 



To obtain contributions from long (C > l e i) trajectories, we require that the dephasing 
length (/^) satisfies > l e i. In this case 7 ~ hence x/\Xl\ scales as R/Ll e i. The 
magnetic susceptibility of two-dimensional ballistic dots (squares) has recently been mea- 
sured by Levy et al. ||. They have obtained x( j^~j°^ ~ kpL ~ 100. In their experiments 
l e i — L — 2L, ~ 3L — 10L . Eq. (|16D , pushed to the limit of its validity, yields [l e i = 1.5L, 
1$ = 8L) x — 250|xl| ; m rough agreement with the experiment. Evidently, a more detailed 
comparison with our theory (e. g., measurements of the crossover field from linear suscep- 
tibility, which is sensitive to disorder) is needed before one may confirm the validity of our 
picture. 
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To summarize, the central point of our analysis was based on the observation that when 
the condition 1$ 3> l e i is met, the weak field behaviour carries the signature of the long, 
disorder dependent trajectories. It is possible that the contribution calculated in Refs. []5]|| 
should be added on top of the present result P^[ . It is the dependence on l e i that cancels 
out which gives rise to the disorder independent result, Eq. (|IfiD, in sharp contrast with, 
e.g., the results expected in the Aharonov-Bohm case. The scale of the field for which the 
susceptiblity decreases may be disorder dependent (Eq. (|17|)) [f26f . 
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FIGURES 



FIG. l.a. A multiply reflected semiclassical trajectory (thick line), and its image within the 
extended zone scheme. A dot denotes the location of a scatterer. 

FIG. l.b. Vector potential within the extended zone scheme (arrow lengths are proportional 
to \A X \). Alternating directions of the staggered field are indicated. 

FIG. 2. The field dependent susceptibility calculated numerically (Eqs. ([!]) and (|TlD ) with 
x = y = for H < Here l d /L = 50.3. Note that MH)) changes sign for H > 

These oscillations will be partially smeared out due to fluctuations in the value of l e \. 
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